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The polarization direction of an electromagnetic field changes and eventually reaches a steady 
state when propagating through a birefringent material with off axis absorption or gain. The 
steady state orientation direction depends on the magnitude of the absorption (gain) and the phase 
retardation rate. The change in the polarization direction is experimentally demonstrated in weakly 
doped (0.05%) Pr®''':Y 2 Si 05 crystals, where the light polarization, if initially aligned along the most 
strongly absorbing principal axis, gradually switch to a much less absorbing polarization state during 
the propagation. This means that the absorption coefficient, a, in birefringent materials generally 
varies with length. This is important for, e.g., laser crystal gain media, highly absorbing and narrow 
band spectral filters and quantum memories. 

PACS numbers: 42., 42.25.Ja, 42.25.Lc, 33.55.+b, 42.55.-f 


INTRODUCTION 

For electromagnetic plane wave propagation in bire¬ 
fringent materials the instantaneous polarization direc¬ 
tion of the electric field vector, E, normally changes dur¬ 
ing the propagation. For example, the polarization of 
a wave propagating along one of the principal axes in a 
(non-absorbing) birefringent material with (initially) lin¬ 
ear polarization at an angle 7 relative one of the other 
principal axes, changes from linear to elliptical (with ma¬ 
jor axes |E| cosy or |E| siny) and then to linear polariza¬ 
tion at angle —y and then back to elliptical etcetera as 
it propagates through the material. For a case where the 
wave is not propagating along one of the principal axes 
the wave vector and the Poynting vector are in general 
not parallel and there are walk off effects (e.g. chapter 6 
in Ref. [l[). 

In this work we consider an initially linearly polarized 
plane wave, propagating along a principal axis in a non¬ 
magnetic off axis absorbing birefringent material. In this 
case polarization rotation can occur even if the initial 
linear polarization is aligned with one of the principal 
axes. In materials with absorption or gain the develop¬ 
ment of the polarization can be complicated. Not only 
could there be oscillations between linear and elliptical 
polarization and walk off effects, there can also be con¬ 
versions from light of one polarization to another polar¬ 
ization due to the absorption (or gain). 

Such effects should for example be present in laser crys¬ 
tal gain media when the gain is anisotropic and the gain 
tensor is not aligned with the principal axes Q. Recently 
there has also been an interest in highly absorbing and 
very narrow bandwidth filters for very specific applica¬ 
tions such as quantum memories Q , or highly absorbing 
narrowband filters with exceptionally large etendue for 
high performance ultrasound optical tomography 30 


and also dynamically tunable high performance filters 0 . 
All these recent papers use filters based on rare-earth-ion- 
doped inorganic crystals where the absorption tensor is 
not aligned with the principal axes. As a consequence 
the light polarization gradually switches during propaga¬ 
tion to a much less absorbing polarization state even if 
the input polarization is perfectly aligned with the most 
strongly absorbing principal axis. This significantly de¬ 
creases the achievable filter attenuation from what might 
be anticipated based on just the material absorption co¬ 
efficient. 

In this work we theoretically analyze the polarization 
rotation in absorbing birefringent materials with a model 
adapted from Q. We provide some simple relations of 
when the effects need to be taken into account and give 
possible suggestions on how to maintain the original high 
absorption. The theoretical analysis is supported by an 
experimental demonstration of the polarization rotation 
effect. 

The paper is organized as follows. First an intuitive 
theoretical background is given. This is followed by 
simulation results and discussions of polarization steady 
states and the effects of absorption on light propagating 
through a crystal. Experimental results are then pre¬ 
sented where the incoming and outgoing polarization di¬ 
rections are studied for crystals with different absorption 
and length. Lastly a remark about the maximum ab¬ 
sorption axis is made before the paper is concluded with 
a summary. A rigorous mathematical treatment of the 
problem can be found in the Supplementary Information 
section. 


THEORETICAL BACKGROUND 

A detailed description of the theoretical framework is 
given in the Supplementary Information, but a simple 
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and intuitive understanding can be achieved by consider¬ 
ing the effects of birefringence and absorption separately. 
In |Fig. 1^ we see the principal axes of a birefringent 
crystal as well as the direction of the transition dipole 
moment, fi. This direction does not coincide with any 
of the principal axes of the crystal. An incoming field, 
Ei„, which is on resonance with the absorbing transi¬ 
tion and with its polarization oriented along one of the 
principal axes, D 2 , generates a polarization, P, along fi 
with a 90° phase shift. This polarization generates a 
field, Ep, along fi with a relative phase relation to Ei„ 
as shown in |Fig. ll o. The total resulting electromagnetic 
field, E(ot = Ei„-|-Ep, now also have a component along 
the Di axis as shown in |Fig. 1| ;. Since the crystal is 
birefringent the polarization state oscillates between lin¬ 
ear dFig. 1| : and |Fig. T} ;) and elliptical ( |Fig. 1^ 1) as the 
field propagates through the crystal. Of course in real¬ 
ity the effects of absorption and birefringence are inter¬ 
twined and occur simultaneously. The theoretical exam¬ 
ination in the Supplementary Information makes it pos¬ 
sible to simulate the propagation of any incoming elec¬ 
tric held polarization through a birefringent crystal with 
an absorption axis with an angle to the principal axes. 
The method works for inhnite incoming plane waves and 
assumes that the lateral (transverse) dimensions of the 
crystal is inhnite as well. Discussions around the results 
from these simulations can be found in the next section. 


SIMULATION RESULTS 


In this section results from simulations of light propa¬ 
gation inside a birefringent crystal with a tilted transition 
dipole moment direction is analyzed. As a specihc ex¬ 
ample a 0.05% Pr^+:Y 2 Si 05 crystal is chosen where the 
transition dipole moment is tilted by 74.6±1.9° from the 
Di crystal axis Q as seen in |Fig. 2] and ca lculated in the 
Supplementary Information, see Eq. (44) 


Input polarization Absorption 
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FIG. 1. (Color online) Simplified view of the propagation of 
an electric field through a birefringent crystal with a transi¬ 
tion dipole moment axis tilted with respect to the crystal axes 
{DiA-D 2 ). (a) The incoming electric field, Ei„, is polarized 
along D 2 ■ (b) The material is polarized and emits an electric 
field Ep along the transition dipole moment which is 180° 
out of phase with the incoming field Ein resulting in the total 
field Etot = Ei„ + Ep which has both a Di and D 2 component 
of the electric field. The magnitude of Ep is exaggerated to 
make the effects easier to see. The electric field is projected 
along the crystal axes Di and D 2 in (c) to make it possible to 
see the effect of propagation through the birefringent crystal 
which causes a relative phase change between the two electric 
field components. In (c) they are oscillating in phase, while in 
(d) they are 90° out of phase and in (e) they are completely 
out of phase (180°). 


riD^ = 1.809 and rib = 1.7851 along the principal axes Q 
is calculated using Eq. (1) to be about 7.25 |xm for light 
propagating along the b axis. 


Phase retardation in a birefringent crystal 

The propagation phase, (j), can be written as (j) = 
radians, i.e. light with the wavelength Aq in vacuum ex¬ 
periences a phase shift of (j) radians when traveling an 
optical length nz through a material, where n is the in¬ 
dex of refraction and z is the propagation depth. The 
difference between the accumulated phases, (j )2 and (j)i, 
for light propagating along two different crystal axes, D 2 
and Di respectively, is termed the phase retardation; 

A , , , „ Anz 

A(/> = ^2 - = 27r—— (1) 

^0 

where An = — njj-^. The crystal length needed to 

obtain a quarter-wave plate, A(j) = for a Y 2 Si 05 crys¬ 
tal with the different refractive indices; udi = 1.7881, 


Polarization steady states 

From the effects explained in |Fig. 1| it is clear that 
it is not possible to maintain a pure linear polarization 
when propagating through a crystal where the transi¬ 
tion dipole moment axis does not coincide with any of 
the optical axes of the crystal. However, in general two 
steady state polarization solutions exist for a forward 
propagating wave. They are elliptically polarized and 
only differ by a 90° rotation of their major axes in the 
D 1 -D 2 plane. The polarization of the steady state does 
not change when propagating through the crystal. The 
ellipticity and the direction of these polarizations depend 
on the following ratio; 
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FIG. 2. (Color online) (a) Transition dipole moment, fi, ver¬ 
sus principal axes for Pr®"'":Y 2 Si 05 (DiTIl 2 T&). (b) Notation 
of linear incoming and elliptical outgoing polarizations are 
seen where 7 is the angle between the Di axis and the incom¬ 
ing polarization and ip is the angle between the Di axis and 
the major axis, a, of the outgoing elliptical polarization. 


R = 


^1)2 ~ ^I?i 

Xahs 


( 2 ) 


where Xabs is the electric susceptibility connected to the 


absorption (for Pr^+:Y 2 Si 05 see Eq. (20) in the Supple¬ 


mentary Information). In general only the steady state 
solution with the lowest absorption or highest gain is sta¬ 
ble. 

In |Fig. 3^ -c the steady state polarization solutions 
when the transition dipole moment makes an angle of 
74.6° relative to the Di axis are shown for three differ¬ 
ent cases where i? ^ 1, i? « 1 and i? ^ 1, respectively. 
For absorption the green curves are stable solutions and 
the red curves are unstable. The reverse is true in the 
case of gain. The first case, seen in |Fig. 3^ , is the same 
as for a 0.05% Pr^+:Y 2 Si 05 crystal which has a transi¬ 
tion dipole moment direction that is tilted 74.6 ± 1.9° 
from Di and absorption coefficients along Di and D 2 of 
aoi = 3.6±0.5 cm“^ and odj = 47±5 cm“^ which cor¬ 
responds to an electric susceptibility for the absorption 
that is Xabs = (8.82±0.8) • 10“^. The difference 
in Y 2 Si 05 is 0.075 which gives a ratio i? ~ 85 1. Here 

one can see that the steady state solutions lie respectively 
along the Di and the D 2 crystal axes, which is always 
the case for i? » 1. 

The next two cases seen in |Fig. 3) 3 and |Fig. '3} : use the 
same values mentioned above except that the absorp¬ 
tion susceptibility, Xabs, is increased by a factor 100 and 
10 000 respectively. In the extreme where i? <C 1, seen 
in |Fig. 3] ;, the solutions lies almost completely along the 
transition dipole moment (unstable for absorption) and 
perpendicular to it (stable for absorption). In between 
these extreme cases when i? Ri 1, seen in |Fig. 3|3 , the 
solutions lie between the crystal axes and the transition 
dipole moments and have a higher ellipticity than in the 
extremes. 

The steady state solutions can be explained by the 


same phenomena that was discussed in |Fig. 1[ i.e. that 
an electric field in D 2 creates an electric field in Di, and 
vice versa, due to the tilted absorption (gain) axis. The 
steady state solutions are polarization states where the 
sum of the components lost and absorbed in Di decays 
at the same rate as the sum of the components created 
and absorbed along D 2 , which in other words means that 
the polarization stays the same. 

The dependence on R can be intuitively explained by 
rewriting the numerator in Eq. (2) as; 


“ ^£>1 = (’^£’2 + nDi){nD2 - noi) 

= ± nDi)^n 


( 3 ) 


If we now assume that udi and are roughly con¬ 
stant then the ratio R measures if An is small/large com¬ 
pared to the absorption and An sets the rate at which 


phase retardation is accumulated (see Eq. (1)). In other 


words if i? ^ 1 then the phase retardation dominates 
over the absorption. This makes all polarizations that 
are symmetric along the crystal axes equivalent (due to 
the effects explained in |Fig. 2] :-e) and therefore the only 
steady state solutions that can exist are along Di and 
D 2 , see |Fig. 3^ . In the other extreme where i? <C 1 the 
absorption dominates and any phase retardation between 
the two electric field components can be ignored which 
leads to one solution orthogonal and one solution parallel 
to the transition dipole moment axis, see|Fig.~3l:. 


Electric field absorption when 1 


From the simulations one can obtain the (expected) 
exponential decay in the electric field components along 
Di and D 2 as a function of propagation. For the exam¬ 
ple of Pr^+:Y 2 Si 05 the decay coefficients can be seen in 
|Fig. 4 1 for a linear input polarization along D 2 (7 = 90° 


Fig. 2 ) 3 ). Even though the input polarization only 


has a D 2 component of the electric field, a Di compo¬ 
nent is almost immediately created in the crystal due to 
the effects explained in |Fig. 1^ -b. The relative size of 
the created component is the same as the ratio of the 
D 2 and Di component of the steady state, since at the 
steady state (when the polarization does not any longer 
change with propagation) the D 2 component is almost 
completely created by the tilted absorption of the Di 
component. Therefore it is a measure of how strong the 
connection between the two components is. For the case 
of pr 3 +:Y 2 Si 05 the initial electric field component along 
Di is approximately a factor of 330 less than the D 2 com¬ 


ponent (see |Fig. 4^ and Eq. (35) in the Supplementary 
Information). 

The exponential decay along D 2 is much larger than 
along Di in the beginning of the propagation due to the 
fact that the transition dipole moment lies much closer 
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FIG. 3. (Color online) Polarization steady states for (a) Pr®''':Y 2 Si 05 with absorption coefficients aoi = 3.6 ± 0.5 cm“^ and 
ao 2 = 47±5 cm~^ which corresponds to an absorption susceptibility Xabs = (8.82 ±0.8) ■ 10“"* and a transition dipole moment 
direction with an angle 8 = 74.6± 1.9° from Di as can be seen in |Fig. 2"^. (b) and (c) has an increased absorption susceptibility, 
Xabs, by a factor of 100 and 10 000 respectively. The ratio R in jEq. (^ is (a) R « 85, (b) R « 0.85 and (c) R « 0.0085. The 
green curves are stable (unstable) and the red curves are unstable (stable) in the case of absorption (gain). 


to I? 2 - However towards the end of the propagation the 
stable steady state solution is reached and both compo¬ 
nents decay at the same rate. Continuing in |Fig. dj j the 
input polarization is changed to 7 = 45° and the polar¬ 
ization steady state is reached sooner since the input po¬ 
larization is more similar to the stable steady state. For 
a linear input polarization the highest absorption rate 
is achieved along D 2 but slows down after propagating 
around 5.4 mm, see |Fig. This leaves us with the ul¬ 
timate limit, which we call the “steady state limit”, of 
the maximum electric field absorption with a high ab¬ 
sorption rate of the incoming light for this material to be 
10“® —10“® and this is reached for a 5.4 mm long crystal. 
To clarify, the light continues to be absorbed after this 
limit but at a much slower rate. In light of this limit it 
is clear that scaling a crystal longer than r; 5.4 mm to 
gain absorption serves limited purpose and one can get 
exponentially higher absorption by avoiding the steady 
state limit. Note also that much earlier, at around 2.7 
mm, the electric field component along Di, which decays 
much slower, is of the same size as the component along 
D 2 and the Di component dominates after 2.7 mm, even 
though the Di component was zero before entering the 
crystal, this component can however be blocked with a 
polarizer after the crystal. 

The (approximate) propagation distance needed to 
reach the steady state limit as a function of the input 
polarization 7 is shown in blue in |Fig. 4| ;. The elec¬ 
tric field component along Di and D 2 when the steady 
state limit is reached is also shown in green and red re¬ 
spectively. Note that the Di electric field component 
is always a factor of ~ 330 larger than the D 2 compo¬ 
nent when the steady state solution is reached. Close to 
7 = 90° the propagation distance decreases quite rapidly 


as 7 departs from 90° and a good control over the pu¬ 
rity of the electric field polarization is required to have 
a high absorption and a steady state limit at 5.4 mm. 
The purity needed is in the order of the ratio of the Di 
and D 2 electric field components for the stable steady 
state, which for Pr®+:Y 2 Si 05 , seen in |Fig. 3^ , is a factor 
~ 330 for the electric field, or equivalently a factor 10® 
for the intensity components (see Eq. (35) in the Sup¬ 
plementary Information). This means that a polarizer 
with better suppression than 10 ® and an angle precision 
of around 0.17° is required to achieve the longest propa¬ 
gation distance with maximum absorption. 

To avoid the steady state limit one has to suppress the 
Di component of the light. This can be done with e.g. a 
linear polarizer. To achieve high absorption it is therefore 
better to use several smaller crystals (lengths of < 5.4 
mm) with linear polarizers blocking any Di component 
in between each set of crystals instead of using one long 
crystal. 


EXPERIMENTAL RESULTS 

To examine the polarization change during propa¬ 
gation several experiments were performed on crys¬ 
tals with different lengths. We used three 0.05% 
Pr 3 -i-:Y 2 Si 05 crystals with different lengths (1 mm, 6 mm 
and 12 mm) and measured the outgoing polarization di¬ 
rection "0 as a function of the incoming polarization angle 
7 , see |Fig. 2 ) 3 . The crystals were kept at cryogenic tem¬ 
perature (2 K). The setup is illustrated in |Fig. 5^ . The 
absorption was probed with attenuated pulses to min¬ 
imize saturation effects on the absorption and between 
each experiment a set of hole eraser pulses were used 
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FIG. 4. (Color online) Shows the decay of the electric field components along crystal axes Di (green) and D 2 (red) as a 
function of propagation distance 2 for a linear input polarization at an angle (a) 7 = 90° and (b) 7 = 45° from Di. The dashed 
black lines are linear fits to the logarithm of the electric field on the form ln(if( 2 )) = ln(i?o) + C ■ z, where the value of C is 
displayed next to each black curve. Note that these values represent the decay of the electric field. To obtain the absorption 
coefficients for the intensity multiply these values with 2, i.e. = 3.6 ± 0.5 cm“^ and aD 2 = 47 ± 5 cm“^. The steady state 
limit can be seen in (a) which shows the largest total absorption that can be reached with a high absorption rate for a linear 
polarization. The inset in (b) is a zoomed in view of the steady state limit reached after « 2.7 mm for 7 = 45°. The oscillating 
pattern has a period of 29 pm (matching a 27r phase evolution calculated with Eq. (1) I and exists due to the interactions 
between the polarization induced in the material and the electric field of the light and the birefringence of the crystal, (c) The 
(approximate) propagation distance needed to reach the steady state limit shown in blue as a function of the input polarization 
angle 7 together with the electric field component along D\ (green) and along D 2 (red) at the steady state limit. Note that 
the a;-axis is cut at 87° and that the axis scale changes afterwards. Note also the plateau reached at 89.83° which shows that 
better angle precision than that is not needed. 


over a wide spectral range to make sure that no perma¬ 
nent hole burning occurred. 

According to the simulations for a Pr^+:Y 2 Si 05 crys¬ 
tal there is a minimum (sufficient) value of the absorption 
where the steady state solution is reached which means 
that the outgoing polarization is almost entirely along 
the Di axis. A 12 mm Pr^+:Y 2 Si 05 crystal is employed 
to illustrate this effect as shown in |Fig. SI d, where three 
measurements at different frequencies within the inho¬ 
mogeneous profile were performed. The inhomogeneous 
line-width, P, of this crystal is measured to be about 9 
GHz. The three sets of measurements at frequencies 0 
GHz, 4 GHz and 13 GHz relative to the line center are 
shown in |Fig. 5| 3 as yellow triangles, cyan squares and 
blue dots respectively together with simulations for each 
case displayed as a red dotted, dashed and solid curve. 
For high and intermediate absorption the outgoing po¬ 
larization is along the Di axis which agrees with the fact 
that the D 2 component is absorbed much quicker than 
the Di component. 

As shown in |Fig. l| ;-e the outgoing polarization shifts 
significantly during propagation. It is therefore impos¬ 
sible to determine if the outgoing polarization is in the 
range ip = 0 ~ 90° or '0 = 90 — 180° for a given input 
polarization without a very precise knowledge of the crys¬ 
tal length. A quarter-wave plate distance was calculated 
previously to be A Ri 7.25 qm which means that a full 
evolution of the phase is done in « 29 pm in Y 2 Si 05 . If 
the crystal length is not known within this precision it 


is indeed impossible to determine the outgoing polariza¬ 
tion direction. The crystal length used in the simulations 
is therefore changed from exactly 12 mm by a few tens 
of pm to the value that best matches the experimental 
results. 

The absorption coefficients in are measured in a 
crystal with a line-width of F = 6.2 GHz and which 
nominally has the same Pr concentration, 0.05% as in 
the present crystal. This would correspond to a peak 
absorption susceptibility of Xal^^ — (8.82 ± 0.8) • 10“^. 
The inhomogeneous profile can be approximated as a 
Lorentzian profile @ that in our experiments have a line- 
width of F = 9 GHz. Assuming that the integral of 
the absorption should be the same for our crystal and 
the crystal used in Q the peak absorption susceptibil¬ 
ity of our crystal should be reduced by a factor 9/6.2 
which gives Xals'^ — (6.08 ± 0.55) • lO”"*. Since the three 
measurements seen in |Fig. 5) 3 are performed at different 
positions on the inhomogeneous profile the absorption 
susceptibility used in the simulations must be adjusted 
accordingly; 

' _ peak (r/2)^ 

Xabs — Xabs (/-/o)'^ + (r/2)^ 

^ Xabs,i = [/ - /o = 0 GHz] (6.08 ± 0.55) • lO'^ 

Xabs,2 = [/ - /o = 4 GHz] (3.39 ± 0.31) • lO'^ 

Xabs ,3 = [/ - /o = 13 GHz] (0.65 ± 0.06) • lO'^ 

( 4 ) 

The simulations agree well with the experimental re¬ 
sults and only in the 13 GHz case is the absorption so 
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low that the outgoing polarization varies from Di. The 
results from the 1 and 6 mm crystals give similar results 
that agree with the simulations. Also here the absorption 
susceptibility must be reduced with respect to the peak 
value to account for the experiments taking place at the 
side of the inhomogeneous profile. 
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MAXIMUM ABSORPTION AXIS 


It has previously been shown, see Ref. [1^, that slow 
light effects split an incoming linear polarization into 
two parts, where one part has a significantly larger slow 
light effect than the other. The transition dipole moment 
direction was then considered as the maximum absorp¬ 
tion direction but it turns out that the maximum (min¬ 
imum) absorption and consequently time delay happens 
when we set the input polarization along the D 2 {Di) 
axis. In general this depends on the ratio R in Eq. (2)[ 
since if i? ^ 1 the maximum absorption direction is al¬ 
most exactly along the transition dipole moment axis, 
but as stated this is not the case for Pr^“'':Y 2 Si 05 where 
R a; 85 » 1 . 



FIG. 5. (Color online) (a) Polarization analyzer setnp. First, 
two polarizers are calibrated against each other by minimiz¬ 
ing the photodiode (PD) signal when there is nothing in be¬ 
tween them. The first polarizer (Poll) and the half-wave plate 
(HWP) makes sure that the incoming linearly polarized light 
is at an angle 7 from Di. The D\ and D 2 directions are mea¬ 
sured with the same setup with the crystal at room temper¬ 
ature where absorption is negligible. By turning the second 
polarizer (Pol2) to maximize the throughput to the photodi¬ 
ode (PD) the outgoing polarization after the crystal is mea¬ 
sured to he tp . (b) Three experiments were performed on a 
12 mm crystal 0, 4 and 13 GHz away from the center of the 
inhomogeneous absorption profile. The simulations have Xabs 
given by Eq. (4) which assumes an inhomogeneous profile in 
the shape of a Lorentzian distribution using a FWHM of 9 
GHz, which was measured in the experiment. The experi¬ 
mental results are displayed as yellow triangles, cyan squares 
and blue dots together with a red dotted, dashed and solid 
curve respectively for the three cases of 0, 4 and 13 GHz. Note 
that the simulation results for the 0 and 4 GHz cases overlap. 


CONCLUSION 

In conclusion, we have modeled and analyzed the po¬ 
larization of light propagating through an absorbing (am¬ 
plifying) birefringent crystal. Steady state solutions for 
the polarization of a propagating wave were found and 
discussed. This led to the conclusion that only increasing 
the length of a birefringent medium (e.g. a rare-earth- 
ion-doped crystal) is not the most efficient method to 
reach higher absorption. This can be an important is¬ 
sue in applications like quantum memories and spectral 
filtering employing birefringent crystals Since the 

model works for any birefringent material with either ab¬ 
sorption or amplification it can be useful in other appli¬ 
cations such as laser crystal gain media or with other 
materials that are not specifically discussed here. 

It also became clear that the polarization direction 
with maximum absorption (gain) depends upon the ratio 
between the phase retardation and the magnitude of the 
absorption (gain). For an absorbing Pr^+:Y 2 Si 05 crys¬ 
tal the maximum absorption axis is along D 2 while the 
stable polarization steady state solution is almost com¬ 
pletely along £> 1 . 

This work will hopefully open up new opportunities 
to investigate more interesting physics and applications 
regarding the propagation effects in rare-earth-ion-doped 
crystals in the future. 
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SUPPLEMENTARY INFORMATION 
Light propagation in a birefringent crystal 

The propagation of an electromagnetic field in a gen¬ 
eral material containing absorption can be a quite com¬ 
plex problem to solve due to the interaction of light and 
matter, but one could reduce this complexity without 
loosing the generality of the solution by considering a 
simplification based on the properties of the medium. We 
start the simplification by considering a plane-stratihed 
medium, which is a valid assumption for the materials in¬ 
vestigated in our case. The assumption of plane-stratihed 
media is a proper approximation for materials that show 
variation in the propagation direction (z-axis) but no 
variation in the lateral {x-y plane) direction. In other 
words, we have homogeneity (same absorption or gain) 
in the x-y plane and possible inhomogeneity (varying ab¬ 
sorption or gain) along the z-axis. The main theoretical 
concept of the calculations in this paper is discussed with 
more details in Ref. 0- 

For macroscopic media. Maxwell’s equations describe 
the dynamics of the helds as follows (the harmonic time 
convention is used): 


I V X E(r,w) = iA:o(coB(r, w)) 

I V X r7oH(r,w) = -j/co(co??oD(r, w)) 


( 5 ) 


where 770 = y ^ is the intrinsic impedance of vacuum, 

Co = 14 /eo/io and fco are the speed and wave number of 
light in vacuum, respectively. In addition, E, B, H, and 
D are the electric held, magnetic hux density, magnetic 
held and the electric hux density, respectively. 

In order to apply Maxwell’s macroscopic equations, it 
is necessary to specify the relations between E, B, H, and 
D. These equations are called the constitutive relations. 


Based on the plane-stratihed assumption, it is natural 
to decompose the electromagnetic held into tangential (x- 
y plane) and normal components. Considering the lateral 
homogeneity and substituting the constitutive relations 
(Eq. (6)) into Maxwell’s equations (Eq. (5) I gives a sys¬ 
tem of ordinary differential equations (ODEs) with the 
variable z. 


dz (ryoJ-Hiyw) 

where J is a two-dimensional rotation dyadic (rotation of 
90° in the x-y plane) as follows: 


J = 


0 -1 
1 0 


and 


J . J = -I 2 = 


-1 0 
0 -1 


( 8 ) 


(9) 


Eq. (7) 


For a more detailed discussion of how to obtain 
which is beyond the scope of this paper, the 
reader is referred to Ref. Q- 

According to the main assumption, the plane-stratihed 
media assumption, the four dyadics e(z), ^(z), ((^( 2 :), and 
/x(z) in Eq. (6) depend only on one spatial variable, 
the propagation distance. The 
which is our medium, is appro- 


z, which describes 
Pr 3 +:Y 2 Si 05 crystal, 
priately modeled as a non-magnetic anisotropic medium 
with zero crossed magneto-electric dyadics. Therefore, as 
a second assumption, we assume non-magnetic proper¬ 
ties (/X = I) and also no crossed magneto-electric dyadics 
(^ = (^ = 0) in Eq. (6) There exists a general set of 


equations for the connection between the fundamental 
matrix M(z) and the four main dyadics in Ref. Q, but 
based on the three assumptions described above and a 
normal incident we can simplify M(z) = M for our case 
as: 


I D(r, w) = eo[e(z, w) • E(r, w) -H r]o^{z, w) • H(r, w)] 
\B(r,a;) = ^[<( 2 , w) • E(r, w) -h ?7o/x(z, w) • H(r,a;)] 

( 6 ) 

where r = xx -I- ?/y -I- zz and w is the angular frequency, e 
and /X are the permittivity dyadic and permeability dyadic 
of the medium, respectively, while ^ and ^ called 
crossed magneto-electric dyadics. All four dyadics and 
all four helds can depend on the angular frequency w, but 
this variable is suppressed below to simplify the notation. 

In a single crystal, the physical and mechanical proper¬ 
ties can often be orientation dependent. When the prop¬ 
erties of a material vary with orientation, the material is 
said to be anisotropic. Alternatively, when the properties 
of a material are the same in all directions, the material 
is said to be isotropic. Bi-anisotropic is a general class 
of linear media which exhibit so-called magnetoelectric 
coupling between the electric and magnetic helds Q. 


where 


M = 


( . “ - - 

\ €.zz 



( 10 ) 





Following Eq. (7)[ we need to investigate the eigenval¬ 
ues and eigenvectors of the M tensor (Eq. (10)) in order 
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to calculate the tangential components of the electric and 
magnetic field. The eigenvectors and correspond to 
the eigenvalues +Ai and — {i = 1,2) for the forward 
and backward waves, respectively: 

for i=l,2 

VO^iy: V0^ix\ 

( 12 ) 

For a full analytical derivation of the eigenvalues and 
eigenvectors see section “Theoretical analysis of the 
propagation matrix M”. 


E-field evolution in the propagation direction 


To understand the E-field and H-field properties while 
the wave is propagating in the z direction, we can em¬ 
ploy the fundamental equation for one-dimensional wave 


propagation (Eq. (7) I. By considering the eigenvectors 
and eigenvalues of the M matrix (Eq. (12)), the wave 
properties along the propagation direction (z) is calcu¬ 
lated for the forward wave: 


dz 


Ui = ikoXiUi Ui(z) = 


(13) 


It is possible to calculate the backward wave, Vi(z), in 
the same way. 

By convention, the polarization of light is described by 
specifying the orientation of the wave’s electric field at a 
point in space over one period of the oscillation. There¬ 
fore, to understand the polarization direction inside the 
medium it is enough to only derive the E-field. The ini¬ 
tial light polarization before interaction with the medium 
could be described by the E-field: 


Now we can extend our fields to a different z using 


Eq. (13) Eq. (15) and Eq. (17) as follows: 


u(z) = -f 


(18) 


Note that since the M matrix does not depend on z 
in this case the eigenvalues and eigenvectors -l-Ai, — Aj, 
Ui and Vi remains constant and it is therefore possible 


to propagate the fields (or equivalently u(z) in Eq. (18)) 


to any desired length in one step. In the more general 
case where the eigenvalues and eigenvectors changes with 
z one has to recalculate these in the simulations every 
distance Az to be able to propagate the fields to z, where 
Az is chosen for numerical stability and resolution. 

Note also that the two polarization steady states for 
a forward propagating wave discussed in the main text 
is simply Ui and U 2 where the solution with the low¬ 
est imaginary part of its eigenvalue is the stable solution 
whilst the other is the unstable solution. This can be 
understood since a high absorption leads to a fast decay 
from that state, i.e. the highest absorption gives the un¬ 
stable state. In the same way, a large negative imaginary 
part corresponding to a high gain will generate a large 
component and therefore become stable since the other 
component can be neglected. 


Light propagation in a Pr®''':Y 2 Si 05 crystal 

We start by calculating the permittivity tensor for 
Pr^+:Y 2 Si 05 and then employ the theoretical approach 
discussed above to obtain the light polarization direction 
while propagating through this specific medium. 

To define the permittivity tensor for Pr^+:Y 2 Si 05 one 
could start by deriving the relation between permittivity 
and susceptibility for the host material (Y 2 Si 05 ) and the 
absorber (Pr) separately as follows: 


Eq — T ^Oyy 


(14) 


In addition, the transverse electric field {x-y components) 
inside the medium could be written as a sum of the E- 
field of the two eigenvectors: 

^inside = AlEi(z) + BE2{z) 

= A{EixX + Eiyy) + B{E2 x^ + E2yy) (15) 

Based on the boundary condition for the tangential 
E-fields at the interface: 


Therefore 


X : Eqx — AEix + BE2x 

y '■ Boy = AEiy -)- BE2y 


A _ EQxE2y—EoyE2x 

Eix E2y — Ely E2 x 

B = — AExx] 


(16) 


(17) 


e = (1 + x) = (1 + Xhost + iXabs) 


(19) 


The imaginary part of the susceptibility is proportional 
to the absorption, while the real part of the susceptibility 
is proportional to the real refractive index. 

In section “Theoretical analysis of the propaga¬ 
tion matrix M” a full analytical derivation of the elec¬ 
tric susceptibility Xabs and the transition dipole moment 


angle 9 is given but the following equations (Eq. (20) 
and Eq. (21)) are approximately correct for the case of 


Pr3+:Y2Si05. 

Xabs is a dimensionless quantity which represents the 
imaginary part of the electric susceptibility and is ap¬ 
proximately proportional to the absorption coefficients 
(a_Di + aD 2 ), in the medium as follows [llj |: 


'^bg^O / , \ 

Xabs = - [O^Di + OiD 2 ) 


LO 


( 20 ) 


where riyg 1.8 is a background refractive index and 

w « 27r-494-10^2 
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The anisotropic absorption coefficients for the two 
sites of a nominally 0.05% Pr^+:Y 2 Si 05 crystal at the 
— ^D 2 transition is measured as shown in the Table 
7.6 in Ref. Q. Based on this measurement the absorp¬ 
tion for site I (605.977 nm) is about 3.6 ± 0.5 along the 
Ill, 47±5 along the ZI 2 , and < 0.1 along the h axis, all in 
the unit cm“^ (see |Fig. "2^ , where T>i, D 2 and b are crys¬ 
tal principal axes). To estimate the susceptibility and 


the transition dipole moment angle 0, Eq. (20) and the 
following equation are used (for more correct values see 
Eq. (44)D; 


%an( 0 )^ 


Qg 2 

CtDj 

= 74.5 ± 1.9° 


( 21 ) 


Xabs = (8.8 ±1.0)-10- 


where Xabs is given in the direction of the transition 
dipole moment shown in |Fig. 2^ . It can be rotated about 
the b axis to the principal axis coordinate system to be 
able to perform the summation in Eq. (19)1 


f Xxx ± iXabs cos^ 0 
iXabs sin 0 cos 0 
0 


iXabs sin 0 cos 0 
1 + Xyy + ^Xabs sin^ 0 
0 


1 


'3.1973 ± 6.2e-5 
2.3e-H 
0 


i 2.3e 
3.2726 ±8.2e-4z 
0 



0 

0 

+ Xzzj 

(25) 


(26) 


Using this permittivity tensor together with Eq. (10)] 
Eq. (12) Eq. (17) and Eq. (18) it is possible to propa¬ 
gate light through an absorbing Pr^+:Y 2 Si 05 crystal. 

THEORETICAL ANALYSIS OF THE 
PROPAGATION MATRIX M 

For light propagating parallel to one crystal axis of a 
plane-stratified, non-magnetic birefringent crystal with 
an absorption (or gain) axis in the transverse x-y plane 
and no crossed magneto-electrical properties the propa¬ 
gation matrix M can be written as; 


crystal _ tj 
Xabs ~ 


-1 . .^dipole 
^abs 


Rz 


cos^ 0 sin 0 cos 0 0^ 
= Xabs I sin 0 cos 0 sin^ 0 0 

0 0 0 , 


( 22 ) 


The next step is to calculate Xhost based on the Sell- 
meier dispersion and the measured coefficients from Ta¬ 
ble 7.1 and Eq. 7.2 in Ref. Q. 


'2.1973 

Xhost I b 
0 



(23) 


Adding up all calculations, the final results for the per¬ 
mittivity tensor in the principal axis coordinate will be: 


(24) 




J 


Aj = ± 


+ iXabs Xabs 


M = 


-njj - iXabs COS^ 0 -iXabs cos 0 sin 0 


I 0 0 -1 0 \ 

0 -1 

0 0 

V -*Xabs cos 6» sin 6» - lyabs sin^ 6» 0 0/ 

(27) 

where riD-^ and Eire the refractive indices of the two 
crystal axes in the x-y plane. 0 is the angle from Di to 
the transition dipole moment axis and Xabs is the electric 
susceptibility associated with the absorption. 

The polarization steady states are given by the eigen¬ 
vectors to the propagation matrix M. Four eigenvectors 
and eigenvalues exists and two of the eigenvalues have a 
positive real part which means they are connected to a 
forward traveling wave. These eigenvalues and eigenvec¬ 
tors are most important for the analysis in this paper. 
The stable eigenvector is the one whose eigenvalue has 
the lowest imaginary part (which is true for both pos¬ 
itive imaginary parts in the case of absorption and for 
negative imaginary part in the case of gain). 


Solving for the eigenvalues of M gives four eigenvalues; 


T 


, Xabs 


— 2i(l — 2 sin^ 0) 


''l >2 "Pl 

Xabs 


1/2 


- 1 


(28) 


Since we are only interested in the forward propagating 
waves only the eigenvalues whose real part is positive are 
of interest. These eigenvalues are denoted Ai and A 2 . 

Before we solve for the eigenvectors in the general case 
we have to account for some special cases where 0 = 0° 


or 0 = 90°. In these cases the transition dipole moment 
axis is parallel to one of the crystal axes Di or D 2 which 
of course results in two steady states, one along Di and 
another along 02- Another special case is when the ab¬ 
sorption is zero, which again is a trivial case where the 
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polarizations of the steady states lie along Di and D 2 . 




and the propagation is given by; 


Now that these special cases have been solved we will 
solve for the eigenvectors in the general case where we 
exclude the special cases mentioned above. The forward 
eigenvectors will be called u^; 

I 

The eigenvector equations we want to solve are; 


Ui(z) = 

vl2{z) = U 2 ( 0 ) e *^“'^^^ 


/ 0 

0 

-1 



( Eix \ 


^ Eix \ 

0 

0 

0 

-1 


E- 

— \. 

Eiy 

- iXabs cos^ 9 

-iXabs COS 0 sin 0 

0 

0 


1 

0 


1 

0 

\ -fXahsCOS0sin0 

-'^'d2 ~ '^Xabs sin^ 9 

0 

0 ^ 


^ 'HoHix j 


^ J 


( 29 ) 


(30) 


(31) 


The first two equations provides a relationship between 
the electric and magnetic fields and are omitted here since 
we are only interested in the ratio between the Di and 
D 2 components of the electric field to get the steady state 
solution, i.e. we only have to solve for the first two com¬ 
ponents of the eigenvector. This leaves us with two equa¬ 
tions; 


in any way we set Eiy = 1 (assuming that Eiy ^ 0). Now 
solving for the ratio Ei^/Eiy we get; 

_ ^ix _ fXabs cos ^ sin ^ (oo\ 

“ ^ ^ ^ +^Xabs cos2 0) 

Using the following ratio; 


(n|,^ + iXabs cos^ 0)Eix + (ixabs cos9sm0)Eiy = \}Eix 

(fXa&s cos ^ sin T fXafas 9) Ely — Eiy 

(32) 

Since we have the option to normalize the eigenvector 

I 


R = 


'-02 




Xabs 


(34) 


Eq. (33) (with Xi from Eq. (28)) can be rewritten in 


the following way; 


Eix 2 cos 9 sin 9 

1-2 cos2 9-iRT sin^ 9f + {iR+{l-2 sin^ 9)f 


(35) 


It can be shown that which means that 

the two steady state solutions only differ by a 90° rotation 
in the D 1 -D 2 plane since the two dimensional rotation 
matrix for an angle of 90° is 


J = (: V) 

Analysis when R —5- 0 

We can now solve for the steady state solutions (i.e. 
electric field components of the eigenvectors u^) in the 


limit i? —^ 0, which in other words means that the ab¬ 
sorption is high compared to the phase retardation (see 


Eq. (3)). The eigenvalues becomes; 



(37) 


i.e. one solution with no absorption (or gain) and a re¬ 
fractive index close to the refractive index of the material 
and one solution with high absorption (or gain) and a 
high refractive index. 

The solution for the ratio of the electric field compo¬ 
nents becomes; 























the eigenvectors in the limit when R ^ oo. 
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Ai : 

gi.x ^ 

Ely 

^ — tan(^) 

A 2 : 

£^2 x ^ 

^ 1 

E2y 

tan 8 


(38) 


Since these are real numbers the polarizations are lin¬ 
ear and the directions are for A2 along the transition 
dipole moment and for Ai perpendicular to it. Since 
the stable solution is the eigenvector whose eigenvalue 
has the lowest imaginary part, Ai is stable for absorption 
and A2 is stable for gain. 


Analysis when R ^ 00 


'Ai ; (n|,^ + ixabs cos^ 8) Eia; + {ixabs cos 9 sin 0)Eiy 

« + iXabs COS^ d)Eia; 

Ely Ri 0 

A2 : (iXabs COS 0 sin 0)E2x + {n\,^ + iXabs sin^ 8 )E 2 y 
~ + iXabs sin^ d)E 2 y 

^ E2x ~ 0 

(40) 

From these equations it is clear that the eigenvector 
corresponding to Ai has almost no Eiy component and 
is therefore along Di and the reverse is true for A2. 


In the other limit where i? —>■ 00 the phase retardation 
dominates over the absorption which gives the following 
approximation of the eigenvalues; 


Ai Ri +iXabs cos^O 

A2 ~ + iXabs sin^ 0 


(39) 


Here the stable solution depends on the transition 
dipole moment direction given by 0. If 0 < 45° A2 is 
stable while if 0 > 45° Ai is stable for absorption and the 


reverse is true for gain. Using Eq. (32) we can solve for 


Analysis of Pr®"'':Y2Si05 


In 0.05% Pr^+:Y2Si05 R is quite high so we can use the 
approximations made in the previous section for when 
i? —>■ 00. Continuing from Eq. (39) we can write the 
imaginary part of Ai as; 



Erom these two equations we can solve for the transition dipole moment angle 9 and the electric susceptibility 
associated with the absorption, Xabs, to be; 


arctan 


(2 • (Im(A2))^+n|,^ ) —"dj 


1/4N 


(2 • (Im(Ai))2-|-n^^) “”di / 

Xabs ~ yj (2 • (Im(Ai))2 + njjy - + ij (2 • (Im(A2))2 -b 


(42) 


In Pr^+:Y2Si05 the exponential decays of the intensity 
along Di and D 2 are measured and they are aoi = 3.6 ± 
0.5 cm“^ and 002 = 47 ± 5 cm“^ respectively. Given 
that the eigenvectors calculated in Eq. (40) are along Di 
and D 2 respectively one can, using Eq. (30)] write the 
following equation; 


can calculate the angle 0 from Di to the transition dipole 
moment axis and the absorption susceptibility Xabs for 
a Pr^+:Y2Si05 crystal using Aq = 605.977 nm, noi = 
1.7881 and 71^2 = 1.809. 


(lmiXi) = ^ 

Im(A 2 ) = ^ (43) 

1^0 = t 


0 =74.6 ±1.9° 

Xabs = (8.82 ±0.8)-10-4 


(44) 


where the factor of 1/2 comes from that we now deal 
with electric fields and not intensities, kg is the wave 
number and Ag is the wavelength of the incoming wave An approximation of 9 and 
in vacuum (not to be confused with Ai and A2 which are 


An approximation of 0 and Xabs in Eq. (42) can be made 
in vacuum (not to be confused with Ai and A 2 which are by assuming that Im(Ai) <C n]j. and that noi ~ nua ~ 
the eigenvalues of M). Using Eq. (42) and Eq. (43) one ribg (which is the case for Pr^+:Y2Si05). This results in; 
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6 ft! arctan 


arctan 


' (Im(A2))^n|)^ 

.(Im(Ai))2n|,^ 


1/4N 


/ lm(A2) 

Im(Ai) 


= arctan 


Q!£)2 

aoi 


(45) 


and 


Xabs ~ 2Im(Ai)n_Di + 2lm{\2)nD2 

_ aOinDi anaDa ^ (Q!Di + <^D2)'n‘bg 


‘‘■0 


^0 




(46) 


Eq. (45) and Eq. (46) are only valid in the case where 


-R —>■ oo but allows us to use Eq. (20) and Eq. (21) for 


the case of Pr^+:Y2Si05 which gives values quite close to 


the more correct values obtained in Eq. (44) 


Continuation of analysis when R —> 0 


In the other case where R —>■ 0 we only have one axis 


with absorption (see Eq. (37) I and therefore in Eq. (43) 
we replace ctDa by adipoie that measures the absorption 
along the transition dipole moment axis for A2 whose 
eigenvector points in that direction and set Im(Ai) = 0. 
9 is now given by the angle from the crystal axis to the 
axis of which the absorption is measured. Assuming that 
Xabs ^ ''t-d- the following proportionality between Xabs 
and adipoie can also be calculated using Eq. (37)] 


^ 1 adipoie 
Xabs ~ 7 : I ; 

2 \ ko 


(47) 


This shows that the electric susceptibility Xabs does 
not always scale linearly with the absorption coefficient 

^dipole ■ 




















